home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / program / nrpas13.zip / CHSTWO.DEM < prev    next >
Text File  |  1991-04-29  |  1KB  |  54 lines

  1. PROGRAM d13r9(input,output);
  2. (* driver for routine CHSTWO *)
  3. CONST
  4.    nbins=10;
  5.    npts=2000;
  6. TYPE
  7.    barray = ARRAY [1..nbins] OF real;
  8. VAR
  9.    glinext,glinextp : integer;
  10.    glma : ARRAY [1..55] OF real;
  11.    chsq,df,prob,x : real;
  12.    i,ibin,idum,j : integer;
  13.    bins1,bins2 : barray;
  14.  
  15. (*$I MODFILE.PAS *)
  16. (*$I GAMMLN.PAS *)
  17.  
  18. (*$I GSER.PAS *)
  19.  
  20. (*$I GCF.PAS *)
  21.  
  22. (*$I GAMMQ.PAS *)
  23.  
  24. (*$I RAN3.PAS *)
  25.  
  26. (*$I EXPDEV.PAS *)
  27.  
  28. (*$I CHSTWO.PAS *)
  29.  
  30. BEGIN
  31.    idum := -18;
  32.    FOR j := 1 to nbins DO BEGIN
  33.       bins1[j] := 0.0;
  34.       bins2[j] := 0.0
  35.    END;
  36.    FOR i := 1 to npts DO BEGIN
  37.       x := expdev(idum);
  38.       ibin := trunc(x*nbins/3.0)+1;
  39.       IF (ibin <= nbins) THEN bins1[ibin] := bins1[ibin]+1.0;
  40.       x := expdev(idum);
  41.       ibin := trunc(x*nbins/3.0)+1;
  42.       IF (ibin <= nbins) THEN bins2[ibin] := bins2[ibin]+1.0
  43.    END;
  44.    chstwo(bins1,bins2,nbins,-1,df,chsq,prob);
  45.    writeln;
  46.    writeln('dataset 1':15,'dataset 2':15);
  47.    FOR i := 1 to nbins DO BEGIN
  48.       writeln(bins1[i]:13:2,bins2[i]:15:2)
  49.    END;
  50.    writeln;
  51.    writeln('chi-squared:':18,chsq:12:4);
  52.    writeln('probability:':18,prob:12:4)
  53. END.
  54.